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Description 

The invention relates to a method of determining the stroke volume and the cardiac output of the 
human heart from the pulse-type blood-stream pressure signal derived from the aorta and consisting in 

5 each case of a systolic and diastolic period. It is Icnown in practice to calculate the stroke volume and the 
cardiac output with the aid of a thermodilution determination by injecting cold liquid into the blood stream 
and measuring it downstream at regular intervals. 

In another known method, a pulse contour method or prescription is used to determine the stroke 
volume Vs and the cardiac output or heart minute volume Q from a pressure signal p(t) measured in the 

10 human aorta. In this connection, the stroke volume is the volume of blood ejected by the heart in one 
contraction or stroke. A typical value is, for example, 70 cm^, but this quantity can vary from stroke to 
stroke. The cardiac output is the volume of blood which the heart pumps in a unit time of one minute. In 
this connection, a typical value is 5 l/min, and this volume can typically vary between 2 and 30 l/min. 

In the past relatively simple relationships have been assumed between the measured aorta pressure 

75 signal p(t) and the mean aorta flow g(t). In this connection, the starting point is the pressure signal since 
said pressure can be measured relatively easily and well, but the flow cannot. In practice it has been found 
that such a relationship or "model" is much more complicated than was first assumed. This has emerged, 
In particular, in hospitals, where this so-called pulse contour calculation has resulted in frustration in the 
case of seriously ill patients. 

20 In practice, usually a so-called "Windkessel" or air-receiver model is assumed, the aorta being 
conceived as a single compliance (figure lb). In this case a windkessel - a container partially filled with 
liquid and with gas - which can absorb the surges in flow and buffer them in the gas bubble and which is 
incorporated downstream of a pulsatile pump is meant. The volume of blood ejected by the heart pump is 
largely received in the aorta or windkessel and partially flows away through the peripheral vascular regions 

25 of the various organs which branch off from the aorta. In the period in which the heart pump does not eject 
a volume of blood, i.e. the diastolic period T^, the outward flow from the windkessel continues to be fed to 
the peripheral vascular regions. The aorta does not then contain any gas but the aorta wall is elastic and 
this elastic vascular wall fulfils the same function as the windkessel. 

This concept contains various factors which are neglected. Firstly, the aorta is much longer than wide. 

30 In fact, the pressure wave generated by the heart requires a time of 0.1 to 0.3 sec. to reach the end of the 
aorta and then approximately the same time to return to the heart, whereas the entire expulsion period of 
the heart (the systolic period or systole. Tg) lasts only 0.2 to 0.4 sec. At the beginning of the output flow, 
only a small portion of the compliance is therefore available for buffering the outward flow. However, still 
more important is the fact that, after two times the propagation or transit time of the pressure wave, twice 

35 the capacity is available, albeit with a doubling of the initial pressure wave amplitude. In addition, a second 
neglect which is generally made is that the windkessel is linear, that is to say a doubling of the stored 
stroke volume is accompanied by a doubling of the pressure. In reality it has been known' for a number of 
years from the work of G.J. Langewouters et al. In J. Biomechanics 17. pages 425-435 (1984) that the 
relationship between pressure p and volume V Is strongly nonlinear In accordance with an arctangent 

40 function (figure 3a). In addition, said arctangent function is dependent on the age and sex of the patient. 
Finally, the drainage from the aorta to the peripheral vascular zones is not concentrated at one location but 
distributed along the aorta. That is to say. it takes a little time before the increased pressure at the start of 
the aorta has reached the various branching points to the peripheral vascular zones and the drainage 
actually increases due to the increased aorta pressure. 

46 To summarise, the neglects therefore amount to the fact that the aorta does not behave linearly and 
that there is a time factor due to a travelling wave along the aorta. 

In practice, however, another model is also assumed, the so-called transmission line model (figure la). 
In this approximation, the aorta is conceived as a homogeneous elastic tube of "adequate" length filled with 
liquid (blood), with drainage to the peripheral vascular zones concentrated at the end. Such a tube has two 

50 characteristic properties. A pressure wave generated by the heart at the start of the tube travels at finite 
velocity Vp to the end of the tube, is partially reflected there and travels back to the start of the tube. By the 
time the reflected wave arrives there, the expulsion phase, the systole, of the heart is already over, the aorta 
valve is already closed and the heart no longer has any trouble from the reflected wave since the latter is, 
after all, held back by the closed valve. The "adequate length" is that length of tube which, given the 

55 propagation velocity of the wave, ensures that the reflected pressure wave does not return too soon. The 
propagation velocity Vp is primarily determined by the area A of the cross section of the aorta, the 
compliance C* of the aorta per unit of length and the density p of the blood in the tube as follows: Vp = 
VA/pC. 
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In addition to the propagation velocity, the second characterising property is the characteristic 
Impedance Zq. This Innpedance reflects the ratio between the anriplitude of the pressure wave and that of 
the accompanying flow wave in the aorta and is also affected by the cross sectional area A and the 
compliance C*. In fact, the characteristic impedance is given by the formula: 

5 

The fact that the ratio between pressure p(t) and flow q(t) is fixed and given by the formula: q(t) = p(t)/Zo 
means that the pressure wave and the flow wave at the start of the aorta are identical In shape at least till 

10 the instant that the reflected pressure (and flow) wave has (have) returned. 

This model also contains various approximations which, although they have less serious consequences 
for the stroke volume calculation, make the latter still too Inaccurate under certain circumstances. Firstly, it 
is also true of this model that the aorta does not behave linearly in relation to the relationship between 
pressure p and volume V and therefore, with constant aorta length, as regards the cross section A = V/l. 

75 Since the compliance C per unit of length is the derivative of the area of the cross section with respect to 
pressure, at the prevailing pressure po: C - [dA/dp]po. A will indeed increase with increasing pressure, but 
at the same time C will decrease because the aorta becomes increasingly less distensible with increasing 
pressure in accordance with the arctangent function mentioned (figure 3b). In a particular pressure region, 
these two factors will compensate for each other in the formula for Zo and the latter will be virtually 

20 constant. At increasingly higher pressures, the decrease in C takes place more quickly than the increase in 
A, and the characteristic Impedance Zo will start to slowly increase in accordance with the root of the 
inverse product of A and C*. The nonlinearity therefore plays a role but this is less serious than in the 
windkessel model. A second approximation is that the aorta is assumed to be homogeneous. In reality, the 
aorta is neither homogeneous in cross section nor in distensibility of the vascular wall. The aorta cross 

25 section decreases towards the periphery. This is initially compensated for because the arterial branches of 
the aorta start to contribute to the total cross section as soon as the pressure wave and flow wave reach the 
branches, in addition, the vascular wall is, however, increasingly less distensible towards the periphery. This 
manifests itself in a decrease in compliance C. Since both A and C in the formula for the characteristic 
impedance decrease over the aorta towards the periphery, Zo will increase. This means that there is no 

30 question of a single reflection at the peripheral resistance at the end of the aorta tube, as assumed in the 
model, but that there are distributed reflections along the entire aorta, initially small but considerably 
increasing later in the systole. With a short systolic period, the heart only experiences the (low) impedance 
at the start of the aorta while contracting, but with a longer-lasting systole it also experiences the (higher) 
peripheral impedance. This will mean that the instantaneously expelled quantity of blood will become 

35 increasingly smaller as the systole advances since it is gradually curbed by the Increasing reflected 
impedance. 

A third approximation is the following. After the blood has been forced into the aorta due to an increase 
in the pressure, the pressure in the subsequent diastolic period will, of course, not be zero immediately, but 
the aorta will slowly empty through the peripheral resistance. As a consequence of this, the subsequent 
40 stroke will be forced out against a somewhat increased pressure and the two pressures will in fact be 
superimposed on each other. The diastolic pressure will always increase further with every stroke until 
equilibrium has been reached between systolic inward flow and diastolic outward flow. This is comparable 
to a windkessel function. 

Finally, it is also true that the drainage is not concentrated at the end of the aorta but is distributed over 
45 the length thereof. 

To summarise, the neglects are a consequence of the fact that the aorta has nonlinear properties and 
cannot be conceived as homogeneous. 

Both the above model concepts are not adequate if an attempt is made to calculate the aorta flow 
pattern from the aorta pressure as a function of time during a systole using the models. The curve found 
50 with the transmission line model only resembles an actual flow curve to any extent under the conditions of a 
short-duration systole. The windkessel model does not actually even permit a calculation of a flow curve. In 
practice, for both models only the stroke volume of the heart, that Is to say the flow curve integrated over 
one systole, is therefore calculated from the pressure curve. 

The object of the invention is to eliminate the abovementioned problems and to provide a method for 
55 the said determination of the stroke volume and cardiac output which takes account both of the windkessel 
properties and the transmission-line properties. This is achieved in the method as defined in the claims. 

In this case, the nonlinear aorta properties are important and are involved in this method, while it is 
advantageous for the verification of the correctness of the steps of the method if the flow curve can be 
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continuously calculated and displayed during a systole. The integral of the flow curve then yields the stroke 
volume. 

The invention will be explained in greater detail on the basis of an embodiment, reference being made 
to the drawings, wherein: 

5 figures la and lb show the two known pulse contour models, viz. the windkessel model and the uniform 

transmission line model; 

figures 2a and 2b show the pulse contour model according to the invention based on the transmission 
line model supplemented by the windkessel function, respectively without and with peripheral drainage; 
figure 3a shows a graph with stylised pressure/cross section and pressure/compliance per unit length 
10 relationship for a 40 year old male aorta; 

figure 3b shows a diagrammatic representation of the calculation of a linearised pressure using the 
arterial final diastolic pressure as initial value; 

figures 4a and 4b show, respectively, an Input pressure signal p(t) and a calculated flow curve q(t) 
calculated with the aid of a pulse contour model; and 
76 figure 5 shows a number of graphs of the pulse contour (line) and thermodllution (o) cardiac output of six 
open-heart operation patients plotted against the number of the thermodilution determination. 
The invention assumes the transmission-line model, in particular, the transmission-line property that 
pressure wave and flow wave are related by a fixed factor, the characteristic impedance Zq, at least before 
reflections return to the heart from the periphery. This model lacks the gradual curbing of the outward flow 
20 from the heart by early reflections and the diastolic windkessel action. According to the invention, these two 
shortcomings are eliminated by means of adding one component, an aorta windkessel compliance Cw, to 
the transmission-line model (see figure 2a). After all, in the systole, the flow is determined in the first 
instance by and is then gradually curbed during the ejection phase by the pressure increase of py,(t) 
which is built up in the windkessel receiver by the inflow q(t). 
26 A further neglect was that the peripheral drainage was conceived as concentrated at the end of the 
transmission line. As a result of now associating, according to the invention, the equivalent peripheral 
resistance Rp with the windkessel compliance Cw, the latter is connected to an average (albeit low-pass- 
filtered) pressure and consequently represents reality better. See figure 2b. For in this model, the outward 
flow 

30 

Qr (t) 
p 

is subtracted from the inflow q(t) so that the curbing of the inflow takes place somewhat less rapidly. 

36 The improved model according to the Invention contains linear timelnvariant components. In reality, we 
are concerned with a strong nonlinear aorta whose pressure/cross section relationship varies in accordance 
with an arctangent. According to the invention two versions are now possible. In the first version, the 
pressure curve recorded in the method is first linearised and then applied to a linear pulse contour model. 
In the second version, the pressure measured in the method is regarded as given, but the pulse contour 

40 model is provided with pressuredependent components. 

The linearised pressure curve model: 

In this first version it is assumed that the pressure curve arises because a certain volume of blood Is 
45 pumped by the heart into the aorta, and this Is followed by a rise in pressure. The magnitude of the rise in 
pressure depends on the pressure/volume relationship of the aorta. At a high pressure, a low volume 
increase will result in a high pressure increase, whereas at a low pressure, the same volume increase 
results in a smaller pressure increase (arctangent relationship). However, It is possible to calculate, by 
means of the pressure/volume relationship assumed to be known, the pressure curve which would have 
50 existed If the aorta had had a linear pressureA^olume relationship. The procedure can be indicated 
graphically in a simple manner (see figure 3b). For every stroke, the start of the expulsion phase in the 
pressure curve and the associated final diastolic pressure Pd are determined. At said pressure, the nonlinear 
aorta pressure/volume relationship or, assuming a fixed aorta length, pressure/cross section relationship is 
linearised by means of the tangent to the curve and the equivalent linear pressure is calculated for every 
55 pressure sample taken on the pressure curve. Assuming an initial pressure pd of 70 mmHg. in this example 
of a pressure sample of 100 mmHg, the equivalent linear pressure becomes equal to approximately 92 
mmHg. This pressure is then applied to the pulse contour model which, for characteristic impedance and 
windkessel compliance calculated from the applicable pressure/cross section relationship, has the values 
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associated with the final diastolic pressure for the stroke concerned. 
The nonlinear model: 

5 In this second version, it is assumed that the pressure curve arises as a result of the outward flow of the 
blood into a nonlinear aorta, given all Its properties and not only its windkessel properties, as above. A 
linear correction is not then applied to the pressure curve but is accepted as it is. However, during the 
expulsion phase of the stroke concerned, the value of the characteristic impedance and of the windkessel 
compliance varies according to the pressure value such as it exists at any instant and is calculated from the 

70 nonlinear arctangent pressure/cross section relationship according to the above formulae. At any instant 
when a pressure sample is taken, the associated characteristic impedance and windkessel compliance, and 
the effect of these altered values on the flow are calculated. 

The use of the arctangent relationship: 

75 

According to Langewouters et al., the relationship between the area of the aorta cross section and the 
aorta pressure can be specified using the formula: 

A(p) = AJ0.5 + (1/^.)tan-i{(p-po)/pi}] (1) 

20 

where A(p) is the cross section as a function of the pressure p, A^ Is the limit cross section at very high 
pressure, and po and pi are, respectively, the pressure at the point of inflection of the relationship and the 
slope of the relationship at the point of inflection. For a typical curve see figure 3a. 

The compliance C* per unit of length as a function of the pressure follows from this as a result of 
25 differentiation with respect to p as: 

C'(p) = [dA/dp]p^ = C„/{l*(p-Po)VPi"} (2) 

30 where 

Cm = A„/^pj * the maximum compliance at p = p^. 

35 The values of the parameters A^, Po and pi are known from the work of Langewouters et al. In particular, 
these parameters depend on the sex and age of the person. The value for Am has, in addition, a dispersion 
of approximately 20%. It is not therefore possible to determine the values of the components in the pulse 
contour model in an absolute sense without more exact calibration unless a standard deviation of 
approximately 20% is regarded as acceptable. 

40 Zo and Vp can readily be calculated with the aid of A and C found in this way. To calculate the total 
windkessel compliance of the aorta, still one more parameter is required, namely the effective length of the 
aorta. The windkessel compliance is then the compliance per unit of length multiplied by the effective 
length: Cw = C A good value for for adults is 60 cm. 
The value of the peripheral resistance. 

45 The pulse contour model of figure 2b includes a component representing the total peripheral resistance 
Rp of all the vascular drainage beds in parallel. The value of the latter is initially unknown. After the stroke 
volume Vs has been calculated for the first time and the heartbeat rate /. the number of heartbeats per 
minute, has been found from the sum of the systolic and diastolic time periods T = Ts + T^,, if / = 60/T. 
the cardiac output can be calculated as the product of the heartbeat rate and the stroke volume for the 

50 stroke concerned: 

Q = A/s. (3) 

The average aorta pressure pm is also known as the integral of the pressure over a time period divided by 
55 the time period: 
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p„ -(i/T) jp{t)dt m) 

T 

5 The value of the peripheral resistance then readily follows as: 

Rp = Pm/Q (5) 

Since the peripheral resistance only varies slowly, its value, calculated for the present stroke, can be used 
70 for the flow calculation of the subsequent stroke. For the first stroke, a reasonable initial estimate will have 
to be taken. Since the value of the flow in the systole is not strongly dependent on the value of the 
peripheral resistance, a rapid convergence is obtained in this way. 
Calculation of the flow using the pulse contour model. 

The flow q at any instant can be calculated using a simulation. In this connection, equations are set up 
76 which describe the behaviour of the model of figure 2b. The aorta flow follows from the equation: q(t) = 
l/Zo).{(p(t) - Pw(t)}. The flow through the peripheral resistance follows from: qR(t) = Pw/Rp; and the volume 
in the air-receiver compliance Cw follows from: Vw(t) = + / {(q(t) - qR(t)}dt. where Vq is the initial 
volume In Cw- Finally, the pressure in the air-receiver follows from the volume/pressure relationship Pw(t) = 
^(Vw(t)). As stated earlier, the stroke volume Vs can be found from the flow curve by integrating the flow 
20 during the systole, the heartbeat rate f as the quotient of 60 and the time period of every stroke, the cardiac 
output Q as the product of stroke volume Vs and heartbeat rate f, and the peripheral resistance Rp as the 
low-pass-filtered quotient of the average pressure pm and the cardiac output Q as indicated in formulae 3, 4 
and 5. 

A measured pressure curve and the flow curve calculated from it are shown in figure 4 for a 62 year old 
25 male patient. Before every stroke, the beginning of the rising edge of the pressure curve in figure 4a is 
determined using a pattern recognition algorithm. At this instant, the diastolic pressure pa on the curve is 
recorded, and with the latter, the correct values for the characteristic impedance and the windkessel 
compliance are found using a function generator and introduced into the model. The pressure is linearised 
using the arctangent pressure/ cross section relationship and is applied to the model. The correct values for 
30 the parameters which determine the arctangent pressure/cross section relationship follow from regression 
equations according to Langewouters et al. and are as follows for women and men, respectively, having an 
age n: 



35 



for women: 


for men: 


Am = 4.12 

Po = 72 - 0.89n 

Pi = 57 - 0.44n 


Am = 5.62 

Po = 76 - 0.89n 

pi = 57 - 0.44n 



where Am is in cm^, the pressures p is In mmHg and the age n is in years. 

Obviously, these calculations can advantageously be carried out by a digital processor. For this 
purpose, it is necessary to replace the continuous pressure curve by a series of samples with regular 
spacing intervals and to repeat the calculation for every sample taken. The flow curve therefore becomes 
available at these time instants. A preferred sample Interval is 10 ms. A shorter interval of 5 ms is possible 
if a fast computer Is used but only gives a marginal improvement in the precision. The integration of the 
flow curve In order to calculate the stroke volume is replaced by a summation of the flow samples in the 
systolic period. 

The precision of the new method. 

The new method has been evaluated in six patients who were undergoing an open-heart operation. The 
reference value for the cardiac output was measured at some time instants during the operation by 
repeating a thermodilution (TD) determination at least four times after injections of cold liquid distributed 
with precise regularity over the respiratory cycle. A single thermodilution determination would be too 
inaccurate. The particular feature is that the four determinations were Injected so as to be spread regularly 
over the respiratory cycle, as a result of which an averaging of the four values usually means an 
improvement in the precision of the estimate of the actual cardiac output by a factor of four. This 
improvement is only achieved with mechanical insufflation and not, therefore, if the patient breathes 
spontaneously. Spontaneous breathing always takes place during the first determination, i.e. at "s". as 
marked in figure 5, and during all the subsequent determinations, the patient Is anaesthetised and he Is 
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artificially insufflated. 

The cardiac output calculated by the novel pulse contour method was recorded at the same time for 
comparison. The results are shown separately for each of the patients in figure 5. The number of 
simultaneous determinations, the mean error, the error standard deviation (SD) and the error span are 

5 shown for each patient number 8910 to 8912 and 8915 to 8917 in table 1 and averaged for the group. Note 
that the sixth patient has a larger error for three determinations. These were carried out after a period of 
total heart-long bypass while, although the heart was already beating and pumping blood, the pump of the 
heart-lung machine was still assisting at approximately one litre per minute. This flow was in fact measured 
by the thermodilution method but not by the pulse contour method, resulting in this discrepancy. 

TO For comparison, the cardiac output was also calculated by the earlier method according to Wesseling et 
al. from US 3,841,313, which method was known as the most accurate in the literature. The mean error 
standard deviation was in this case 12%, against 8% for the novel method. The novel method therefore has 
a variance of less than half that of the old. 

Roughly speaking, all the mean errors - shown in the last two rows - with the new method are more 

75 than a factor of V2 lower than with the older method. 



Table 1 



20 



25 



pat nr 


TD 


PC 


Error old 


Error new 






old 


new 


mean 


SD 


span 


mean 


SD 


span 


8910 


4.67 


5.21 


4.53 


0.54 


0.43 


1.02 


0.14 


0.38 


1.12 


8911 


5.43 


4.89 


4.99 


0.54 


0.38 


1.02 


0.44 


0.37 


0.99 


8912 


4.59 


5.31 


5.41 


0.73 


0.45 


1.24 


0.08 


0.38 


1.37 


8915 


4.10 


4.49 


3.80 


0.38 


0.65 


1.81 


0.30 


0.26 


0.60 


8916 


4.69 


4.98 


5.13 


0.29 


0.45 


1.44 


0.43 


0.33 


0.89 


8917 


4.29 


4.33 


4.03 


0.04 


0.88 


2.97 


0.26 


0.50 


1.50 


mean 


4.63 


4.87 


4.50 


0.24 


0.54 


1.58 


0.13 


0.37 


1.08 


SD 


0.46 


0.39 


0.52 


0.45 


0.19 


0.74 


0.30 


0.08 


0.33 



During the intensive monitoring of a patient in the operating theatre or Intensive-care ward, it is only 
rarely possible to record the arterial pressure in the aorta. For reasons of convenience and safety for the 
patient, the pressure in the artery radialis or brachialis or a finger pressure will normally be recorded non- 
invasively. This pressure Is distorted compared with the aorta pressure, that Is to say the shape of the 
pressure curve has changed. This is due to the fact that the pressure curve generated in the aorta 
propagates through the arterial system to the more peripherally situated upper arm, pulse or finger and that 
this transmission path has certain properties, the most important of which are its narrow-band nature and its 
low-frequency resonance. It is possible to measure these properties in a large number of patients. In this 
connection, the inventor has found that the transmission properties have a relatively small spread around a 
mean value. In that case it is possible to correct the resonance properties by an anti-resonance In a filter, as 
a result of which it is indeed possible to accurately determine the pulse contour stroke volume of the left- 
hand half of the heart. To a good approximation, a simple second-order filter can be used for this purpose. 
In practice, a digital, IIR (infinite impulse response) second- or third-order filter is used which corrects the 
mean resonance precisely up to a frequency of 8 Hz. 

In this case, on the basis of a series of Input signal samples Xn, a series of output signal samples yn is 
calculated using the formula: 



50 



M 
I 
k= 



<ikyn- 



where c and d are constant filter coefficients. For the second-order filter, the following values are used: 
Co = 0.077 
di = 1.484 
d2 = -0.561 
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Claims 

1. A method for determining blood stroke volume from a pulsatile aorta blood-flow pressure signal in a 
subject, comprising the steps of: 

5 - providing a sensor for measuring arterial blood pressure in the subject and producing a 

corresponding aorta blood-flow pressure signal consisting of a systolic and and a diastolic period, 

- providing a simulation model of the aorta as a transmission line supplemented with a windkessel 
compliance, and providing an arctangent relation for the nonlinear pressure-volume relation in the 
aorta, 

10 - introducing the pressure in above model and calculating the flow therefrom as function of time, 

and 

- integrating the flow over the systolic period thereby providing the stroke volume. 

2. Method in accordance with claim 1, wherein said step of providing the arctangent pressure-volume 
75 relation is derived, with a fixed aorta length, from the aorta pressure/crosssection according to: 

A(p) = M0.5 + (1/ir).tan-i{(p-po}/pi}]. where 
A(p) = the cross section as a function of the pressure p, 
Am = the limit cross section at very high pressure. 

Po = the pressure of the point of inflection of the arctangent relationship, and 
20 Pi = is the slope of the arctangent relationship at the point of inflection. 

3. Method in accordance with claim 2, further including 

- continuously introducing the measured pressure into the simulation model, and 

- continuously adapting the characteristic impedance of the transmission line and windkessel 
25 compliance to the pressure of the windkessel compliance in accordance with the arctangent 

relation during each stroke. 

4. Method according to claim 2, further including 

- linearizing the pressure on the basis of a linear pressure/volume relation of the aorta, and 

30 - introducing said pressure into the pulse contour model while maintaining the values of char- 

acteristic impedance of the transmission line and windkessel compliance associated with the final 
diastolic pressure. 

5. Method in accordance with claim 4, further including calculating the flow g(t) according to a formula 

35 q(t) = (l/Zo).{p(t)-pw(t)}, where = characteristic impedance = Vp/AC), and Pw = the 

compliance pressure, and 
p = density of blood, 
C'(p) = C„^{1 +(p-Po)'/pi'}. where 
Cm = Am/wpi = the maximum compliance at p = Po- 

40 

6. Method in accordance with claim 1, further including calculating cardiac output Q as product of the 
stroke volume and heart rate. 

7. Method in accordance with claim 6, further including 

45 - calculating a peripheral resistance Rp in accordance with a formula Rp = pJQ, where pn» = 

mean pressure = (1/T)./p(t)dt, and 

- calculating the flow through the peripheral resistance Rp in accordance with a formula 



50 



Qr (t) = P„/Rp. 



8. Method in accordance with claim 1 . further including 

- measuring the blood pressure in the arteria radialis or brachialis, and 
55 - filtering the measured pressure signal in order to correct resonance distortion in the pressure 

curve caused by transmission. 

9. Method in accordance with claim 1 , further including 
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- measuring the btood pressure, non-invaslvely, in the finger, and 

- filtering the measured pressure signal in order to correct resonance distortion in the pressure 
curve caused by transmission. 

5 Patentanspriiche 

1. Verfahren zur Bestimmung des Blut-Schlagvolumens aus dem pulsformigen Aorta-Blutstromungs- 
Drucksignal in einem Subjekt, umfassend die Schritte: 

- Vorsehen eines Sensors fOr die Messung des arteriellen Blutdrucks in dem Subjekt und Erzeu- 
10 gung eines entsprechenden Aorta-Blutstromungs-Drucksignals, bestehend aus einer systolischen 

und einer diastolischen Periode, 

- Vorsehen eines Simulationsmodells der Aorta als Ubertragungsteitung, erganzend mit einer 
Windkessel-Compliance, und vorsehen einer Arcustangens-Beziehung fGr die nicht-lineare 
DruckA/olumenbeziehung in der Aorta, 

15 - Einfuhren des Druckes in das obige Modell und Berechnen der Stromung daraus als Funktion der 

Zeit, und 

- Integrieren der Stromung uber die systolische Periode, was das Schlagvolumen ergibt. 

2. Verfahren gemaB Anspruch 1, wobei der Schritt des Vorsehens einer Arcustangens- 
20 DruckA/olumenbeziehung abgeleitet wird, mit einer festen Aortalange, aus der Aorta- 

Druck/Querschnittbeziehung entsprechend: 
A(p) = AJ0.5 + (1M.tan-^{(p-po)/pi}], worin. 

A(p) = Querschnitt als Funktion des Drucks p, 

A„, = der Grenzquerschnitt bel sehr hohem Druck, 
25 po = der Druck an dem Wendepunkt der Arcustangensfunktion, und 

pi = die Steigung der Arcustangensfunktion am Wendepunkt. 

3. Verfahren gemafl Anspruch 2, auBerdem enthaltend 

- kontinuierliches Einfuhren des gemessenen Druckes in das Simulationsmodell. und 

30 - kontinuierliches Einpassen der charakteristischen Impedanz der Obertragungsleitung und der 

Windkessel-Compliance an den Druck der Windkessel-Compliance gemaB der Arcustangersbe- 
ziehung wMhrend jedes Schlags. 

4. Verfahren gemSB Anspruch 2, aufierdem enthaltend 

35 - Lineahsieren des Drucks auf der Basis einer linearen DruckA/olumenbeziehung der Aorta, und 

- EinfOhren des Drucks in das Pulskonturmodell, wahrend die Werte der charakteristischen Impe- 
danz der Obertragungsleitung und der Windkessel-Compliance. zugeordnet dem finalen diastoli- 
schen Druck, aufrechterhalten werden. 

40 5. Verfahren gemaB Anspruch 4, auBerdem enthaltend das Berechnen der Stromung q(t) entsprechend 
der Formel 

q(t) = (1/Zo).{p(t)-pw(t)}. wobei Zo = charakteristische Impedanz = yJplAC , Pw = der Compliance- 
Druck, und 
p = Blutdichte 
45 C'(P) = C^{1 +(p-po)2/pi2}, wobei 

Cm = Am/wpi = die maximale Compliance bei p = po. 

6. Verfahren gemaB Anspruch 1 , auBerdem enthaltend die Berechnung des Herzminutenvolumens Q als 
Produkt des Schlagvolumens und der Pulszahl. 

50 

7. Verfahren gemaB Anspruch 6, auBerdem enthaltend 

- Berechnen eines peripheren Widerstandes Rp gemSfl einer Formel Rp = p^/Q, wobei Pm = 
mittlerer Druck = (1/T)./p(t)dt, und 

- Berechnen der Stromung durch den peripheren Widerstand Rp entsprechend einer Formel 

55 
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8. Verfahren gemafi Anspruch 1, auSerdem enthaltend 

- Messen des Blutdrucks in der Arteria radialis Oder brachialis, und 

- Filtern des gemessenen Drucksignals zwecks Korrektur einer Resonanzverzerrung in der Druck- 
kurve, die durch Transmission verursacht ist. 

5 

9. Verfahren gemaB Anspruch 1 , aufierdem enthaltend 

- Messen des Blutdrucks non-invasiv im Finger, und 

- Filtern des gemessenen Drucksignals zwecks Korrektur einer Resonanzverzerrung in der Druck- 
kurve, die durch Transmission verursacht Ist. 

10 

Revendications 

1. M^thode pour determiner le volume battement sanguin h partir d'un signal de pression de d^bit 
sanguin aortique puls§ chez un sujet. comprenant ies Stapes consistant h : 
15 - fournir un capteur pour mesurer la pression sanguine arterielle chez le sujet et produire un signal 

de pression de d^bit sanguin aortique correspondant constltu^ d'une p^riode systolique et d'une 
periode diastolique, 

■ fournir un modele de simulation de I'aorte en tant que ligne de transmission, complete par une 
composante compliance de Windkessel, et fournir une relation arc tangente pour la relation non 
20 lin^aire pression-voiume dans I'aorte, 

- introduire la pression dans le module ci-dessus et calculer le d^bit h partir de celul-ci en fonction 
du temps, et 

- integrer le d^bit sur la periode systolique en fournissant ainsi le volume battement. 

25 2. M^thode selon la revendication 1, dans laquelle ladite ^tape consistant k fournir la relation arc tangente 
pression-voiume est d^riv^e. pour une longueur d'aorte fixe, de la pression/section transversale de 
Taorte selon : 

A(p) = A„,[0.5 + (1M.tan-H(P-Po)/pi}],ou 
A(p) - la section transversale en fonction de la pression p, 
30 Am = la limite de la section transversale k pression trds ^lev^e, 

po = la pression du point d'inflexion de la relation arc tangente, et 
Pi " est la pente de la relation arc tangente au point d'inflexion. 

3. M6thode selon la revendication 2, comprenant en outre Ies Stapes consistant ^ 
35 - Introduire contlnuellement la pression mesur^e dans le modele de simulation, et 

- adapter continuellement Timpedance caracteristique de la ligne de transmission et la compliance 
de Windkessel a la pression la de la compliance de Windkessel selon la relation arc tangente 
durant cheque battement. 

40 4. Methode selon la revendication 2, comprenant en outre Ies etapes consistant h 

' lin^ariser la pression sur la base d'une relation lin^alre pression/volume de Taorte. et 

- introduire ladite pression dans le modele du contour de la pulsation tout en maintenant Ies 
valours d'imp^dance caracteristique de la ligne de transmission et de compliance de Windkessel 
associees h la pression diastolique finale. 

45 

5. Methode selon la revendication 4. comprenant en outre une etape consistant a calculer le debit q(t) 
selon la formule 

q(t) = (l/Zo).{p(t)-Pw(t)}, oO Zo = rimp^dance caracteristique = Vp/AC, et Pw = la pression de 
compliance et 
50 p = densite du sang, 

C'(p) = CJO + {p-Po)2/pi2}.oCi 

Cm = A„/7rpi = la compliance maximale ^ p = po. 

6. Methode selon la revendication 1. comprenant en outre une etape consistant k calculer le debit 
55 cardiaque Q comma produit du volume battement et du rythme cardiaque. 

7. Methode selon la revendication 6, comprenant en outre Ies etapes consistant k 
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- caiculer une resistance p^riph^rique Rp selon une formule Rp = Pn/Q. oD pm = la pression 
moyenne = (1/T)./p(t)clt, et 

- caiculer le d^bit via la resistance p^riphdrlque Rp selon une formule 



qRpW PvN-^p. 



10 8. M^thode selon la revendication 1, comprenant en outre les etapes conslstant h 

- mesurer la pression sanguine dans I'artere radiale ou brachlale* et 

- filtrer le signal de pression mesur^ afin de corriger la distorsion de resonance dans la courbe de 
pression due h la transmission. 

75 9. Methode selon la revendication 1 , comprenant en outre les etapes consistant a 

- mesurer la pression sanguine, de mani^re non invasive, dans le doigt, et 

- filtrer le signal de pression mesur^ afin de corriger la distorsion de resonance dans la courbe de 
pression due k la transmission. 

20 
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